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Abstract 

We examine critically the Gambier equation and show that it is the generic linearisable 
equation containing, as reductions, all the second-order equations which are integrable 
through linearisation. We then introduce the general discrete form of this equation, the 
Gambier mapping, and present conditions for its integrability. Finally, we obtain the 
reductions of the Gambier mapping, identify their integrable forms and compute their 
continuous limits. 



1. Introduction 



The classification of the integrable second order differential equations, based on 
their singularity properties, resulted to four classes [1]: 

- equations that are simple derivatives of integrable first order equations 

- equations that are autonomous (i.e. they do not have any explicit dependence on 
the independent variable) and which are integrable in terms of elliptic functions 

- equations which are non autonomous but in which the the independent variable 
appears in some simple form (linearly or at most quadratically) and which define 
the P transcendents 

and finally 

- equations which are non autonomous and in which the independent variable enters 
through some free functions. These equations are solved by linearisation i.e. they 
can be converted to a linear differential system. 

Prominent among this last class is the Gambier equation [2]. This equation is, in 
fact, the generic equation of the linearisable class, in the sense that all the others can 
be obtained as its special limits. The essential feature of the Gambier equation is that 
it describes the coupling of two Riccati equations in cascade (i.e. the solution of the 
first Riccati equation appears in the coefficients of the second one). Thus the Gambier 
equation is best written as: 

y' = -y 2 + by + c, (1.1a) 
x = ax 2 + nyx + a, (1-1°) 

where a, b and c are functions of the independent variable, a is a constant which can 
be scaled to 1 unless it happens to be and n is an integer. The precise form of the 
coupling is indicated by the singularity analysis which, moreover gives constraints on 
the coefficients a, b and c. In fact, out of these three functions only two (in general) 
are free. Eliminating y between (la) and (lb) one can write the Gambier equation as a 
second order ODE: 
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x = h a xx + bx a x + (a — abjx 
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An important remark is in order at this point. The equations of the Painleve/ Gam- 
bier classification are usually given in canonical form, which means that all possible 
transformations of the dependent and independent variables have been used in order to 
simplify their form. This does not seem to be done in the case of the Gambier equation. 
Indeed, as we will show in the next section, a suitable transformation of the dependent 
and independent variables allows us to put 6 = 0. Thus, the Gambier equation contains 
only two functions, which moreover are constrained by the integrability requirement. 

The discretisation of the Gambier system leads naturally to what we have called the 
Gambier mapping. In [3] we have proposed such a discretisation which we have studied 



using the discrete analog of the singularity analysis, namely the property of singularity 
confinement. In this paper we propose to reexamine the discrete form of the Gambier 
equation and determine its most general expression. Once this form is established we 
can proceed to the study of its particular, limiting, forms and propose expressions for 
the remaining linearisable discrete equations. For the sake of completeness we calculate, 
in the next section, the various limits of the Gambier equation in the continuous case. 

2. The Gambier equation and its various limits 

The canonical list of second order equations with the Painleve property is still an 
unsettled question. The simplest way out of the dilemma is to adopt the attitude of 
Gambier [2] who has presented a minimal list of 24 equations which contain, in principle 
all the basic equations. The remaining ones can be obtained through what in modern 
parlance would be called Miura transformations. Among the equations of the Gambier 
list some belong to the linearisable family. Here they are: 

(G5) x" = -3xx' -x 3 + q(x' + x 2 ), 



(G13) x" = h q q + rxx' + r'x 
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(G14) x" (l -) '/•'••'•' , ' •; (2.1) 
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n + 2 nx 
To this list one must, in principle, add the equation 

(G6) x" = -2xx' + qx' + q'x 

which is nothing but the derivative of the Riccati equation. It is easy to show that the 
Gambier equation (G15) contains all the previous ones: it is in some sense the general 
linearisable equation. Instead of using (G15), which corresponds to a = 1 in (1.2), 
we will work with (1.2) itself where one can directly see the relation to the coupled 
Riccati' s. 

First we start with (1.2) for a = 1, and reduce it to its canonical form. For this we 
introduce the following transformation of the independent variable t to a new variable 
T through ^ = g where g is defined by = an d simultaneously X = gx. 

This leads to an equation of the form (1.2) with 6 = 0, which must be considered its 
canonical form (similarly (G15) is canonical for = 0). Moreover the Painleve property 
requirement introduces one further relation between a and c (or, equivalently, between 
/ and V0- 

Equation (G14) is the easiest to obtain: it suffices to take a = 0. The canonical 
form corresponds to b = c = 0. Indeed, in addition to the independent variable trans- 
formation which allows to put 6 = 0, when a = we have an additional gauge freedom 



which allows us to put c = 0. 

„ n — 1 x' 2 n + 2 , a 2 q , 9 . 

x" = + a xx x 3 + a'x 2 . 2.3 

n x n n 

(with b = c = equation (1.1a) leads to y = z ^ zo and (1.1b) for a = is reduced to a 
linear ODE for 1/x). 

Equation (G13) requires that we take the limit n — > oo on (1.2). The result is 
(where d = lim cn) : 



x' 2 x' 
x" = — + axx' + bx' - —a + (a - ab)x 2 + dx - bo. (2.4) 



Equation (2.4) is (G13) in non canonical form. In order to reduce it to the standard 
expression we take 6 = and introduce a gauge x — > px such that d = p" / p — p' 2 /p 2 . 
The equation reduces then to 



x = h q q + rx x + r x . (2-oJ 

x x 



which is just (G13). What does the limit n — > oo really means in the level of the coupled 
Riccati's? Since n goes to infinity y must go to zero for the equation to remain meaning- 
ful and thus the quadratic term in (1.1a) disappears. The canonical form corresponds 
to b = and a new function is introduced through d = nc. Finally, if we divide (1.1b) 
by x and take the derivative a term ny' appears, which from (1.1a) is equal to d. Thus, 
equation (G13) is, in fact, nothing but a derivative of a Riccati after we have divided 
by the dependent variable. 

Finally in order to obtain (G5) we start by taking n = 1 which makes the x' 2 /x 
term vanish. Integrability implies a = and we choose a = — 1, c = 0. This leads to 
the equation: 

x" = -3xx' - x 3 + b(x' + x 2 ) (2.6) 

which is (G5) in canonical form. Finally it seems that (G6) is not in any sense related 
to the Gambier equation (G15). 

3. The discrete analog of the Gambier equation, revisited 

The discretisation of the Gambier equation is based on the idea of two Riccati 
equations in cascade. The discrete form of the first is simply: 

_ ay + b 

where y = y n and y = y n +i- The denominator of (3.1) can be generically be brought to 
this form through a scaling of y and a division by an over-all factor. The second equation 
which contains the coupling can be discretised in several, not necessarily equivalent, 
ways. In [3] we have proposed the discretisation: 

x = fXV + a (3.2) 
1-gx 



A different approach could be based on the direct discretisation of (2.3) in the form: 



x — x = —fxx + (gx + hx)y + k (3-3) 

In what follows we shall not choose a priori a particular form. We shall rather start (in 
the spirit of [4]) from a generic coupling of the form: 

axxy + (3xx + yxy + 5x + exy + Qx + r/y + 9 = (3-4) 

Implementing a homographic transformation on x and y we can generically bring (3.4) 
under the form: 

xx + yxy — exy — 9 = (3-5) 

(the sign changes were introduced for future convenience). A choice of different transfor- 
mations can bring (3.4) to the form (3.3) while (3.2) can be obtained through a special 
choice of the parameters of (3.4). Note that (3.4) contains an 'additive' type coupling 
xx + Sx + (x + rjy + 9 = for special values of its parameters, but the generic form 
(3.5) is that of a 'multiplicative' coupling where 7, e do not vanish. Solving (3.5) for x 
we obtain the second equation of the discrete Gambier system in the form: 

exy + 9 , 
x = — . 3.6 

x + 7y 

Clearly, a scaling freedom remains in equation (3.6). We can use it in order to bring it 
to the final form: 

x -= x y' d+ / (3.7) 

x + dy 

Eliminating y and y from (3.1), (3.6) and its upshift, we can obtain a 3-point mapping 
for x alone but the analysis is clearer if we deal with both y and x. 

The main tool for the investigation of the integrability of the Gambier mapping 
will be the singularity confinement criterion [5]. A first remark before implementing 
the singularity confinement algorithm is that the singularities of a Riccati mapping are 
automatically confined. Indeed, if we start from x = (ax + 13)/ (jx + 5) and assume that 
at some step x = —5/7, we find that x diverges but x and all subsequent x's are finite. 
Thus, the intrinsic singularities of (3.6) do not play any role. However, the singularities 
due to y (obtained from (3.1)) may cause problems at the level of (3.6). Whenever y 
takes a value that corresponds to either of the two roots ±c of the equation: 

y 2 - c 2 = (3.8) 

we obtain x = ±c/d irrespective of the value of x and thus the variable x loses a degree 
of freedom. On the other hand, once we enter a singularity there is no way to exit 
it unless y assumes again a special value after a certain number of steps. Thus, if we 
enter the singularity through, say y = c we can exit it through y = — c after iV steps. 
However, if y were to take the value c again some steps after taking it for the first time, 



then it would take it periodically and the singularity would be periodic. This is contrary 
to the requirement that the singularity be movable: a periodic singularity (with fixed 
period) is 'fixed' in our terminology. 

The first singularity condition can thus be obtained in the following way. We 
assume that at some step y assumes the value c solution of the condition (3.8). This 
value of yo = cq evolves under the action of the Riccati and we obtain after iV steps, 
yN ■ We require that 

Vn = -cat, (3.9) 

i.e. the second root of (3.8). It is thus straightforward to write the first confinement 
conditions for the first few values of N. We have for instance 

ac + b 

N = 1 + c = 

N = 2 ,( a c + g) + + 6(c + l) + - = < 3 ' 10 > 
ac + o + c + 1 

and so on. The equivalent of this requirement in the continuous case is that the res- 
onance be integer. We see here that the discrete condition is much more complicated 
and while one can easily compute the first few instances no general expression can be 
given. Once y passes through the second special value — c, there is a possibility for x 
to recover its lost degree of freedom through an indeterminate form 0/0. This is the 
confinement condition. In full generality (and somewhat abstract form) it reads: 

xn + d N y N = (3.11) 

or, using (3.9), 

xn = d^CN (3-12) 
where a; at is the iV-th iterate of x through (3.6). We have for example 

N = 1 - d = cd 

AT o 1 (ac + b)c + ddc 2 (c + l) = = (3-13) 

Jy = 2 = = = dc. 

d c(c+l) + dd(ac + b) 

The two confinement conditions put constraints on the coefficients a, 6, c and d just as 
the Painleve requirement restricts the parameters in the continuous case. The better 
approach is to start with given a, c and use (3.10) to solve for b. The second condition 
becomes then an equation for d. For iV = 1 we find explicitely dd = c/c. For N = 2 
the equation for b is linear and the one for dd is just a homographic mapping with 
coefficients depending on a, 6, c. 

One important question that remains to be addressed is that of the continuous 
limit of the Gambier mapping. We start from the system 

v = j^t (3 - 18a) 

x -= X y' d+ / (3.186) 
x + dy 



and introduce the following expansions for the parameters 



a = 1 + (a H )e, 

7 n 

^ = 7e 2 , (3.19) 

C= T £ ' 

d = 1 + Se, 



and for the dependent variables 

" = ra £ (3 ' 20a) 

where f = 5 + 7'/ (27) and ft = /'//. We obtain at the limit e — > the two Riccati's: 



X' = -f 2 X 2 + nXY + 1 (3.21a) 

aft ft 2 ft' 

— + — -- 

n n z n 



where the coefficient of the coupling term n = 2c /b is a priori a function but with hind- 
sight we have ignored its derivatives. While the continuous limit takes quite expectedly 
the form of two Riccati's in cascade we have still to show that they are indeed of the 
Gambier form and in particular that the coefficient of the coupling term n is in fact an 
integer and equal to N. 

The key to this proof is the first confinement condition. Let us start with yo = c. 
Given the dependence of c on e (3.19) we have 

Vo = ^e 2 (3.22) 

In order to do away with the e 2 factor we introduce the auxiliary quantity i/j through 
y = e 2 ip and rewrite (3.22) as: 

^0 = ^ (3.23) 



The confinement condition is 



= (3.24) 



let us now compute ifjjsr using the discrete Riccati (3.18) at lowest order in e. Substituting 
the expressions (3.19) of a, b we have at lowest order: 

^ = ^ + 7 (3.25) 

Thus ipN = ip + N*y and substituting the values of ip and ipN we find n = N. Thus 
the coupling coefficient does indeed go over to the integer iV which is the number of 
steps required for confinement. (Had we started from yo = —c we would have obtained 
n = —N. The fact that ±c play different roles is due to the fact that the discrete Riccati 
(3.18) is not symmetric with respect to the upward-downward evolution). 



We must remark here that the above continuous limit is incompatible with N = 1. 
Indeed for N = 1 condition (3.13) implies 2^5 + 7' = which would make (3.20) 
meaningless. This is related to the fact that in the continuous case n = 1 is never 
integrable for a = 1. In fact for N = 1 the only meaningful limit of (3.18) is the one 
that takes a to zero. It turns out that this limit is different from (3.20). In fact taking 

V = ye (3.26a) 

and a = 1 + (a + — )e, with 6, c, d as in (3.19) with n = 1, we obtain at the limit e — > 
the system: 



A" = XY (3.27a) 
y' = -Y 2 - aY + 7 (3.276) 

We remark that in this case the equation for X becomes almost trivial since the X 2 
term vanishes together with a. 

4. NONGENERIC FORMS OF THE GAMBIER MAPPING 

An exhaustive study of all the nongeneric cases of the Gambier mapping is a task 
that lies beyond the scope of this work. In principle one has to go back to the system 
(3.1)-(3.4) and, following the steps of the derivation of (3.5), identify all instances where 
some transformation cannot be applied. The bulk of the resulting calculations makes 
this problem hardly tractable and we prefer, in what follows, to limit somewhat our 
scope. 

We start thus with the Gambier mapping in its reduced form (3.1)-(3.7) and 
consider the cases where the coefficients that have been assumed to be nonvanishing, 
do vanish. We are thus led to the systems given below. The equation for x assumes one 
of the following forms: 

(«) 

+ (4.2) 

X 

x= X -^-. (4.3) 



while that for y is given by: 



_ ay + b 

V = (4.4) 

(4.5) 

y 

y = y + b (4.6) 



all the other cases obtained from (3.1)-(3.7) can be brought to one of the above using 
homographic transformations on x and y. Next we shall investigate the singularity 
confinement property of the system consisting of one of the (4.1), (4.2), (4.3) coupled 
to one of the (4.4), (4.5), (4.6). There exist in principle 9 possible couplings, the one of 
(4.1) and (4.4) being the full discrete Gambier system studied in the previous section. 

In order to investigate the coupling (4.1)-(4.5) we apply the singularity confinement 
method. The principle is the same as for the full Gambier case: we enter a singularity 
when y passes through the value c. In order to confine this singularity we require that, 
after N steps, y pass through — c and moreover x assume an indeterminate form 0/0. 
The condition for y to be equal to — c can be worked out for the first few values of N: 

N = 2 8 ' K '" + ' = 4 
ac + 1 

The corresponding conditions for the denominator of x to vanish (which, in view of 
(4.7), entails the vanishing of the numerator) read: 

N = 1 2 = ^ 

lc(ac + l) + ddc 2 c _= ( 4 - 8 ) 

N = 2 =— J =cd. 

d c 2 + dd(ac + 1) 

In order to obtain the continuous limit of this system we start with the equation for 
y. The only continuous limit of (4.5) is obtained for y = i + eY. However this is 
incompatible with the integrability condition where y assumes the values c and — c 
(after N steps). Thus the system (4.1-5) although integrable as a discrete system does 
not possess an integrable continuous limit. 

Next we consider the coupling (4.1)-(4.6) which can be treated just as the previous 
case. The first few conditions for the singularity to be confined are: 

N=l c + b=-c,_ 

N = 2 c + b + b=-c. y > 



combined with 



N = 1 -=cd 



d 

lc(c + b)+c 2 dd _= ( 4 - 10 ) 

N = 2 =-i J = cd. 

d c + dd(c + b) 

In this case the continuous limit is obtained through: x = eX, c = €7, d = 1 + eS and y = 
Y . From the constraint (4.9) on 6, c we find that at lowest order we have b = —2c/N. The 
continuous limit is then straightforward, but one must also verify the second integrability 
condition (4.10). It turns out that the resulting form is noncanonical. In order to bring it 
under canonical form a further transformation is needed on X: X = 7(1 — 5W) / (1+5W) 
and, moreover, we must take 7' = 0. In the case N = 1, (where we have from (4.10): 
26 + 7V7 = 0), the canonical form can be recovered through a simpler transformation, 



X = -7 + 1/W, leading to the linear equation: W" + W'Y/'j + W(Y/^)' = which 
for Y = reduces to just W" = 0. 

We turn now to the case of the mapping (4.2) coupled to any of the three no- 
mographic for y (4.4-6). A general remark is in order here. The mapping (4.2) has as 
only singularity y = oo, i.e. x is defined independently of the value of x only when 
y = oo. Once y in mappings (4.4) or (4.5) hits this special value, (4.2) loses one degree 
of freedom and cannot recover it because y cannot become infinite again (unless the 
mapping for y is periodic which we have excluded from the outset). Thus the combi- 
nation of (4.2) with either of (4.4) or (4.5) is never integrable. On the contrary (4.2) 
coupled to (4.6) is always integrable because the latter, being linear, can never lead to 
y = oo. In this case we find at the continuous limit equation (G6). Indeed, putting: 
x = 1 + eX, y = 2 + e 2 Y and b = /3e 3 , c 2 = -1 + 7c 2 we obtain X' = -X 2 + Y + 7, 
Y' = (3. Eliminating Y leads to X" = —2XX' + (3 + 7' which can be brought to the 
canonical form (G6) through a simple translation. 

Analogous arguments do apply to the case of the mapping (4.3). The singularity 
of this mapping occurs only if y = 0. Again, the argument of y taking twice the value 
being possible only if (4.4-6) are periodic, precludes the integrability of (4.3) coupled to 
any of these three. However there exists a case where y cannot vanish. This is the case 
of (4.4) for 6 = 0. This is the only integrable case of mapping (4.3) coupled to (4.4). 
However it is a trivial one. By transforming y — > 1/y both mappings become linear. 

5. Conclusion 

In this paper we have examined the Gambier equation in both its continuous and 
discrete forms. For the continuous Gambier system we have shown that it is the generic 
second-order differential linearizable system: the other second order linearizable ODE's 
can be obtained as special limits of the Gambier equation. In the discrete case we 
have obtained the Gambier mapping starting from the most general discrete Riccati in 
cascade system (instead of introducing an ad hoc parametrisation as we did in [3]). This 
most general form has made possible the interpretation of the number of steps necessary 
for confinement. In the particular case of the Gambier mapping this integer coincides 
with the one appearing in the coupling term of the ODE's obtained in the continuous 
limit and which is equal to the resonance of the Painleve expansion. 

This remark raises two important issues. The first is whether there exists a system- 
atic relation between the Painleve resonance and the number of steps for confinement, 
i.e. the length of the singularity pattern. We believe that the answer is, in general, 
negative, despite some tempting results like the Gambier system. The second remark 
is even more crucial. Since the Gambier system confines for any number of steps N, 
the limit N — > 00 does in principle exist. We can thus wonder what is the meaning 
of confinement that requires an infinite number of steps. How can one distinguish the 
N = 00 confining case from a nonconfining one? Although we cannot offer a rigorous 
statement, we can present some elements of an answer based on our experience with 
integrable discrete systems. In a nonconfining system the analysis of the singularity 



shows that there is no possibility for confinement ever. In many cases one can even 
formulate this in rigorous terms and prove the impossibility of confinement. In the case 
of a confining mapping the analysis indicates that the possibility of confinement does 
exist but is simply delayed (and pushed to infinity at the limit). More complicated situ- 
ations may exist, those, among others, involving the discrete derivatives of homographic 
mappings. Clearly, at this level the refinement, the notion of confinement itself becomes 
quite delicate. 

The reduced cases of the discrete Gambier system have been only cursorily studied 
in this work. The particular case where the Gambier mapping reduces, for N = 1, to 
the 2-dimensional projective system was not contained in the forms studied here. In 
order to obtain it one must go back to the initial complete form of the discrete Gambier 
system and perform the appropriate reductions there. This question is under active 
investigation [6]. 
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